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A post-Minkowskii approximation of general relativity is described as a power series expansion 
in G, Newton's gravitational constant. Material sources are hidden behind boundaries, and only the 
vacuum Einstein equations are considered. An iterative procedure is outlined which, in one complete 
step, takes any approximate solution of the Einstein equations and produces a new approximation 
which has the error decreased by a factor of G. Each step in the procedure consists of three parts: first 
the equations of motion are used to update the trajectories of the boundaries; then the field equations 
are solved using a retarded Green function for Minkowskii space; finally a gauge transformation is 
performed which makes the geometry well behaved at future null infinity. Differences between this 
approach to the Einstein equations and similar ones are that we use a general (non-harmonic) 

, gauge and formulate the procedure in a constructive manner which emphasizes its suitability for 

■ implementation on a computer. 

o\ . 
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I. INTRODUCTION 

^ . 

Close binary systems are a possible source for the LIGO or VIRGO gravitational wave detectors, but the non- 
' hnearities of the Einstein equations have made such systems difficult subjects of analyses. Although, most notably, 
the post-Newtonian expansion [Q-Q has now been used to calculate the gravitational wave form resulting from the 
inspiralling evolution of the binary through terms of order {v/cf 
. In anticipation of the day when such analytical results can be pushed no further, we have developed a purely 
On ' constructive, iterative approach that allows for straightforward numerical implementation. We present a variation of 
, a post-Minkowskii expansion for rclativistic systems. Formulated as an iterative procedure, it takes an approximate 
solution to the Einstein equations and produces a better one — the error decreases by a factor of G with every 
iteration. The procedure reproduces the standard lowest order post-Newtonian equations of motion and the traditional 
O quadrupole formula for gravitational waves. But higher orders are generally too complicated for analytical work 
^ '■ without the slow-motion assumption of the post-Newtonian approximation. However, on a computer the nth iteration 
5-H of the process is no more difficult than the first, and it is there where we envisage putting this procedure to good use. 

Two sources in a close binary system are likely to contain regions of extreme curvature. However, when the sources 
^ , are far enough apart the multipole moments of the individual sources depend very little on the presence of a companion 
and the evolution, say, of two neutron stars differs only slightly from that of two black holes or any other small, massive 
. , objects. In our approach a boundary surrounds each source — and we focus on the vacuum Einstein equations in the 
' region outside of these boundaries. 
■ ■ ' In many ways our approach is a combination of methods developed by others. It starts with a formal expansion 
of the Einstein equations in powers of Newton's gravitational constant, G, in a manner similar to Kerr Q. The first 
order results are formulated using a multipole expansion in terms of symmetric trace free tensors described by Pirani 
, Thorne |^ and Blanchet and Damour [0 . Particularly at first order our results formally appear quite similar to 
Thome's |^ analysis of linearized gravity, but our multipole expansion is done about the moving boundaries which 
hide the sources. At higher orders, we mimic the approach of Blanchet and Damour but allow for a general 

(not harmonic) gauge which is restricted only enough to be well behaved at future null infinity. Thr oug hout our 
analysis, we rely heavily upon the rigorous mathematical results provided by Blanchet and Damour pQ-na. 

One original aspect of our approach is the freedom from gauge requirements — at least within the constraints of the 
metric being considered a tensor field on a flat background space-time. Also original is our analysis of the equations of 
motion which, at lowest order, is similar to that of Bel et al. except that we enclose the sources within boundaries 
and avoid any formally divergent integrals. At each order of the approximation we find the equations of motion as 
consequences of the desire to have the world line of the center of the boundaries actually match the trajectory of the 
physical source which is contained inside. 

Also, we often use flat-space outgoing-null spherical coordinates attached to the accelerating world lines. These 
provide globally well behaved coordinate systems and seem particularly well suited for problems involving black holes 
as long as the relevant accelerations are small, mv <C 1. 
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In §11 a summary of our assumptions and approximations is presented and followed by descriptions of the notation, 
coordinate systems and formulation of the general Einstein equations which we use. §111 gives a description of the 
iterative procedure along with specific details of how to start the process at the first order, and how to iterate the 
field equations at nth order. But, for the field equations to have a well behaved solution it is first necessary to iterate 
the equations of motion as described in The behavior of the gravitational field at future-null infinity is discussed 
in §0. Some details are relegated to the Appendix including a description of the retarded Poincare transformation — a 
convenient generalization of the Lorentz transformation which relates outgoing-null spherical coordinates centered on 
different world lines. 

Our process is described in a manner that should make implementation on a computer straightforward, particularly 
for the analysis of the binary inspiral problem from the time of the post-Newtonian applicability perhaps down to the 
innermost stable circular orbit, where speeds could be ~ c/2. We expect this approach to fail when the evolutionary 
time scale is comparable to the dynamical time scale. 



II. BASIC FORMULATION 



A. Assumptions and approximations 

We assume that Newton's gravitational constant, G, is small, but not necessarily very small, when the units are 
such that a characteristic mass and distance are of order 1 for a specific problem. This is essentially equivalent to 
assuming that the metric of space-time deviates only modestly from being flat. We also use units where the speed 
of light is unity — the characteristic time corresponds to the light travel time across the characteristic distance. We 
specifically do not assume that the speeds are small, but we do assume that accelerations are, v°' = 0{G). This is 
consistent with the assumptions of weak fields — when fast, not-strong- field sources interact their accelerations are 
small because the fields are weak. 

We deal only with the vacuum Einstein equations and assume that any material sources are shrouded by boundaries 
at surfaces of constant r in outgoing-null spherical coordinates centered on each source. Some geometrical data can 
be given on the boundary in order to distinguish, say, a system involving black holes from one of neutron stars. This 
distinction may be difficult to implement on a computer, and numerical results may only be of general validity and 
unable to carefully examine features of the geometry which depend upon the detailed nature of the sources. But, at 
least in principle, if this process converges to an accurate solution of the Einstein equations, then the detailed nature 
of the true physical geometry at the boundary can be examined to see whether it is physically consistent with any 
particular source of interest. 

We have only modest gauge restrictions which are imposed to keep the geometry well behaved at future null infinity. 
In particular, we do not require the harmonic gauge. 



In §IV we assume that the mass monopole moment. A, is larger than all others by 0{G). This is not a required 
assumption, but it is reasonable on physical grounds and allow us to describe the dynamical equations in familiar 
terms. 

And in §0 we assume that there exists an initial spacelike hypersurface in the past of which the geometry is an 
exact solution of the Einstein equations. The geometrical data on this hypersurface must satisfy minimal gauge 
requirements at large values of r. And the first order approximation to the Einstein equations, h^^, must match up 
smoothly with the data on this initial hypersurface. 



B. Notation 



This formalism considers a gravitational field as a tensor field on a flat Minkowskii background. The mathematical 
notation is that of flat-space tensor analysis with Minkowskii coordinates {t, x, y, z), along with the usual flat metric, 
rjab, and its inverse, 7]°"^, which are both (—1, 1, 1, 1) down the diagonal and zero elsewhere. The operator Va is just 
the usual derivative operator of flat space, and = VaV". Only rjab and ?7°^ are used to lower and raise tensor 
indices, and a tensor index always refers to a component in Minkowskii coordinates. 

For a binary system the positions of the two sources are described by two different world lines of the form z°'{s) 
parameterized by the Minkowskii proper time, s. In the discussion the focus is usually upon only one world line at a 
time with four- velocity, v"". It is convenient to use = z°'{s). 

Occasionally we use outgoing-null spherical coordinates {s,r,9,(j)) based on the world line z'^{s). The scalar field 
s is defined at an event, V, of Minkowskii space by s('P) = s(Q), where Q is the vertex of the future null cone from 
z'^(s) which contains V; it is convenient to use Sj, = s(a;"). In a similar manner v"", as well as any other tensor defined 
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along the world line, can be promoted to a tensor field over all space-time by parallel transport of v°'{Q) up the future 
null cone of Q to T'. 

Any field derived from a tensor defined along the world line, has a particularly simple expression for its derivative. 
For example 

Vbv^ - -hv" (1) 

where 

ka = -VaS, (2) 

and a dot denotes a derivative with respect to the retarded time s, so that 

v" = v^Vbv'' = dv^/ds. (3) 

The quantity r, in the outgoing-null coordinates, is the Minkowskii spatial distance between Q and V as measured 
in the instantaneous rest frame of the world line at Q, 

Ka;-^) ^-t'aK- 2-^(2)]. (4) 
Pirani |^ shows that is the null vector field 

fc'^ = {x'' - z%Q))/r, (5) 



pointing from Q to V . Also 



and 



k'^Va = -1, (6) 



= rib + rUkaij" (7) 
where rib = kj, — Vb is an outward-pointing, spatial unit vector. It is useful to know that 

rVah = rjab + Vakb + Vbka - kakb{l + rkcv"). (8) 

The angles 9 and (/) at V are defined in the usual way with the origin at Q and with a set of orthonormal basis 
vectors which is Fermi- Walker transported along the world line and parallel transported up the null cone. 

The projection operator onto the spatial three manifold instantaneously orthogonal to v'^ at Q is f'^'' = rj"-^ + v°'v^, 
and the alternating tensor orthogonal to is e"^"^ = e°-^^'''Vd- But, note that when is promoted to a tensor field 
via parallel transport up the future null cone it does not become the spatial three metric of a constant t surface if the 
world line is accelerating. 

The description of tensor multipole moments often requires a large set of, say, I indices; we follow Blanchet and 
Damour [ p^ and define a tensor multi- index, ^ = di...di to denote a succession of I space-time indices. We use 
= n'^^ . . . vf-^ for the tensor outer product of / vectors, n'^, and Vl = V^i ■ . • V^, for a succession of / derivatives. 
Often we sum over I from to oo, and this summation is assumed to converge without justification being given. 
Sometimes a set of tensor indices are symmetric, spatial with respect to w"(s) and completely trace free, these are 
referred to as being SSTF. If a tensor has all of its indices SSTF and is a function of only s, then it is denoted by a 
capital, script base letter. Damour and Iyer [Q give a host of useful formulae for decomposing a tensor into SSTF 
parts. We follow their notation and equivalently denote the SSTF part of a tensor by = A"^^^ = ^<<'i - ai>^ 
Also, [1/2] is just the largest integer less than or equal to 1/2. 

C. The Einstein equations on a flat background manifold. 

A metric, g^b, on a four dimensional space-time may be considered as a two indexed, symmetric invertible tensor 
field on a flat, background Minkowskii space. It is convenient to define /i"'' by 

-gab ^ ^ab _ f^ab ^g) 
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and an Einstein tensor density as a functional of 

E^^h) = i-g)i2R''^ - g'^^R), (10) 

so that the vacuum Einstein equation is 

E'^^ih) = 0. (11) 

Landau and Lifshitz ||l^ give an exact form for we write this as 

E^^ih) = -V^/i"^ + VVc/i'^'' + V^Vc/i'^'' 

- Tj''''VcVdh'"^ - 167rr°''(/i) (12) 

where 

167rT"^(/i) = -2W ch"^"^ dh!^" - h^^Vc^dh'"^ 
+ h^'^Vd^ch"'' + h'^^Vd^ch"" + Wdh^'^Wch'"' 

- /i^-^VeVrf/i"^ + Vch^'^Vdh'"' + 16ni-g)Tll{h). (13) 

The quantity T^^{h) is the Landau-Lifshitz pseudotensor, Eq. (96.9) in ||l^ or Eq. (20.22) in Misner et al. ]l6[ |. 
The Bianchi identity translated onto the flat background takes the form 

V.i?'"' = (r;,Tl - ^L)i?"^ (14) 

where F^^ is the usual Christoffel symbol for the real space-time metric, gab, and is 0{G). 

III. ITERATIVE PROCEDURE 

We formally expand the gravitational field in powers of G: is zero, and at first order hf = O(G^) and 
exactly matches the geometrical data on the initial hypersurface. We iteratively assume that i?"''(/i„_i) — 0(0"^), 
with no gauge restrictions on (in particular it need not be in the harmonic gauge), and look for a correction, 

Shf = Ci(G"), such that 

h-\x)^h-'L,{x)+SK\x;G) (15) 
and E"-''{hn) — 0(0"+^). The dependence of 6h'^ on G is allowed to be more complicated than just being proportional 



to G", and in Eq. (15) that functional dependence is explicit — usually the dependence on G is just understood. 
For a given h'^_^^ the next order approximation follows from a solution of 

V^J/^f = i;"^(/i„_i) + 0(G"+i), (16) 

for 5h'^ = 0{G"), with the additional restriction that 

Va(5/if = 0(G"+i). (17) 

That /i-J*!^! + 6h'^ is a more accurate solution to the Einstein equations is revealed by substitution into Eq. ( |T^ ) 
resulting in 

E^^hn-l + 6h„) 
- [E'^'iK^, + 5K) - E''\K-i)] + V^SK" + 0(G"+i) 

-16^K^(/i„_i + SK) - T-\K-i)] + 0(G"+i). (18) 

The first equality follows from Eq. (|l^) , and the second from Eq. ( p^ . With the restriction (^^ each of the first three 
terms on the right hand side are 0(G"+^), and T°'^{h) is quadratic in /i"'' and its derivatives, so the fourth term is 
also 0(G"+i). Thus 

E''\K-i+6K) = 0{G'^+^), (19) 
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and one full step of the iteration consists of solving Eq. ( |16[ ) with the restriction (|17|) . 

The restriction should not be considered a gauge condition. After all, at the nth step /iJ^'Li need not satisfy 
any particular gauge choice. And the restriction does not involve the residual value of ^ ah'^_i in any manner. Thus, 
at nth order there is no limitation upon Vah^''. Also, the divergence of Eq. ( p^ with the Bianchi identity ( p^ implies 
that V^VqJ/iJ^'' = ©(G""*"^); thus, with the proper choices of initial data and of boundary conditions the restriction 
follows naturally from the wave equation (^6|). To emphasize finally the gauge freedom allowed, note that after each 
iteration is complete a A"j = 0(6*") gauge transformation can change the metric by hf^'' h'^^ + dX^ where 

aAf EE V^A,^, + V''A^-77'^^VcA^. (20) 

This transformation preserves the accuracy of the approximation and changes E'^^{h„) only at 0{G'^'^^) and only 
through the change in T'^^{hn). Also at the nth order, the metric may be changed by the addition of a small, arbitrary 
symmetric tensor, ^'^j^i = 0(G"+^); h"^ —> h"^ + 7^+1 only changes E°-^{hn) at 0(G"+^). Such changes are used in 
§0 to insure proper asymptotic behavior. 

As an iterative procedure this is slightly more general than the post-Minkowskii expansion of Blanchet and Damour 
[0. If any initial approximation to the Einstein equations has E°-^{hi) — 0{e) for some small epsilon, then after one 
step E°-''{h2) — 0{eG). And if h'^'' represents an exact solution, independent of gauge, then the procedure terminates. 

The remainder of this paper focuses on a specific, constructive method for performing one full step. 



A. First order approximation 



We formally require at first order that hf' must match the geometrical data on the initial hypersurface. However 
this limitation is only used in in addition, for applications on a computer we are unlikely to use exact initial data. 
For these reasons in this section we just look for an ft,"'' which resembles two moving sources and has E°'^{hi) ~ 0{G^) 
but does not necessarily match onto good initial data. 

A general multipole source, Mi''^{s), confined to a world line, z°'{s), has an ft"'' ~ Sh'^'' which both satisfies 

00 „ 

V^ftf = -47r^ / Mf^{s)VLS'^ix' Zs)ds (21) 
;=o '' 

and also is of the general form (see Appendix ^ 

00 „ 

hf^Yl / M^^^LGix - z,) ds 

00 

= ^VL[r-iMf^(s)]. (22) 
/=o 

We write Mf''-^ as a sum of terms involving completely SSTF tensors, which are 0(G) and functions only of s, in a 
manner which parallels Thome's Eq. (8.4) or Blanchet and Damour's Eq. (2.25) [0. The resulting most general 
multipole source has 

The reverse parentheses imply symmetrization on a and b. Note that of these SSTF tensors, 6f , Cf , and 
always have / > 1, while and Gi have I > 2. 

But the ftf from Eq. (||) can represent the first term in the expansion of the gravitational field of the multipolar 
source of Eq. (l2^ ) only if the restriction (|l^) is also satisfied. A lengthy analysis of the consequences of restriction (|lj) 
invokes a useful identity described in Appendix]^, Eq. (B2), liberally uses /""^ = rj'^'^ + v°'v'^ and V^G = —4ttS'^{x — Zg) 
and ultimately results in 

oc ^ 

Vaftf (x) ^Y. [^''-^^ + '^1^ + ^''q'^'Cf'^ + f'^'Vi'^ + v'^Bi + f'"^'£^-^ + Tl^ + (\'^'gf~^ 

1=0 •' 

_^fbcU^L-l ^ ^fcpL ^ ffbL _^ ^b^d,jqL-l ^ fd, j^L-l^g ^Ql^^ _ _^ 0{G^). (24) 
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The completeness of the decomposition of SSTF vectors and tensors allows the decomposed parts in Eq. ( p^ to be 
matched up according to the location of the index — whether it sits on v'', e^'^'^' or on an SSTF object. Thus 
the restriction ( p7| ) requires that the multipole source, Mf**^, satisfy 

+ =0{G^) [vX (25) 

+ + Tif = 0{G^) [SSTF''], (26) 

+ + Jfi^ - 0(G2) (27) 

+£f +A:f =0(G2) [/'"^']. (28) 

These are closely related to Eqs. (8.5) of Thorne [§ and Eqs. (2.26)-(2.28) of Blanchet and Damour 0. 

The multipole moments Ai . . . /Cf represent the distribution of a source as computed with respect to the world 
line And our Ai . . ./Cf are closely related to similar quantities used by Thorne The differences are that our 
Bi, Ci and I?f correspond to the negative of his; our Q^, Hi and Ji correspond to twice his; and our multipoles 
are defined with respect to an arbitrary world line, his are with respect to the origin of the Cartesian coordinates. 
But this close correspondence easily allows us to follow his analysis for physical interpretations of some of the lowest 
multipole moments. In particular (^i + I)i)w°/4 + (S° + 7ii)/4 corresponds to the four momentum of the source 

relative to the world line; ^['^(^J' + v'^ — Hi)/2 + e'^q''{C'[ + corresponds to the total angular momentum of 

the source about the world line; —{Ai -\-'Df — l-Li)/{Ai + 2?i) corresponds to the displacement of the center of mass, 
away from the world hne; and (C" + ji)/2 corresponds to the spin angular momentum. Also and G{ each give 
a transverse (to n° and w** at large distances) trace-free contribution to hf' and, therefore, represent the moments of 
sources of gravitational waves which are inside the boundaries. These interpretations reflect the algebraic resemblance 
of these multipole moments, calculated via on a boundary close to a source as it moves along an accelerating 
world line, with the corresponding physical quantities as usually defined at large distances in linearized gravity. 

Thorne describes a gauge transformation which can be used to set all of I?f , 7i f , and /Cf to zero while still 
preserving both Eqs. (^6|) and ([l7|). In this gauge, then, any choice of the ylf , , Cf , £f , and which satisfies 
Eqs. (|25|)-(p8|) could represent the first order multipolar decomposition of a gravitational source. A simple, interesting 
source of this type has only Ai and non-vanishing to give the source both mass and angular momentum but no 
additional structure. 



B. The nth order approximation 



In this section we discuss the procedure by which the next order approximation is found, but we avoid issues of the 
behavior of the approximation at large r — that subject is analyzed in § [v| . 

A given {n — l)th order approximation to the Einstein equations has E^_i = E°-'^{hn-i) = 0{G'^). A step in the 
iterative procedure requires a solution of Eq. ([l^), with 5h'^ = 0{G") which satisfies the restriction (|l^). Following 
Blanchet and Damour ||l^,|ll|], we represent Sh'^'' as the sum of a particular solution of the inhomogeneous wave 
equation, pf^, and a general solution of the homogeneous wave equation, q'^^; thus, 

Sht^pt + qt, (29) 

where 

v^K' - K'-i (30) 

and 

V'^Qn'' = 0, except on the world hne. (31) 

And q'^'' is chosen so that 

VaK' + Va'zf = 0(G"+i). (32) 
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1. The solution for p"^^ 



The quantity p"^ is formally given, with a retarded Green function, by 

= / K-i{^')G{x - x') d^x'. (33) 

In a numerical implementation this evaluation is the single computationally intensive element. It would most likely 
be performed through standard finite differencing of the inhomogeneous wave equation. The fields would propagate 
freely through the interiors of the inner boundaries where the source, E^_i, would be set to zero. Boundary conditions 
for p'^ would only be imposed in the wave zone. It is convenient to interpret p'^ as resulting from the nonlinearity of 
the Einstein equations with /iJ^'Li creating an effective stress-energy outside of the boundaries. Similarly, 5^'' appears 
to come from an nth order correction to the multipole moments of the source hidden behind the boundary. However, 
with a nonlinear theory, the split of dh!^ into two parts in a manner which depends upon the location of the boundary 
is at least modestly arbitrary, and these interpretations are only suggestive. 



2. The solution for q'^^ 

As a solution of the homogeneous wave equation, except on the world line, q^^ has a general representation similar 
to that of hi'' given in § |[II A| , 



qt{x)^Y.^L[r-'M!^'\s,)]. (34) 
1=0 

The majority of this section shows how the SSTF components of Af are determined by the SSTF components of 
E!l^_i on the inner boundaries; the results are exhibited in Eqs. (^7|)~(^0|). 

A consequence of the formal solution for p"^ in terms of the retarded Green function is that 

VaK" = ^ / V'jKf_i(x')G(x - x')] d^x' 

/ Vl(i?f_i)G(x-x')dV, (35) 

where use is made of integration by parts and the symmetry of the retarded Green function. From the Bianchi identity 
(H) and F^^ = 0(G), the second integrand is Cl(G"+i), and 

VaPf = ^ 1 V,[E^'.,{x')G{x - x')] d^x' 

+ 0(G"+i). (36) 

This four- volume integral reduces to boundary integrals about each of the two sources and a third at large r', which 
gives a vanishing contribution, for x fixed, as r' goes to infinity — this can be verified by a lengthy analysis which 
starts with the multipolar expansion of the Green function provided by Blanchet and Damour ||lffl. 

From this last equation it is clear that outside of the inner boundaries the 0(G") part of VaP^ is a homogeneous 
solution of the vector wave equation. And we must obtain its multipolar decomposition in order to find the corre- 
sponding (7°^, a homogeneous solution of the spin- two wave equation which also happens to satisfy Eq. (^2|). The 
needed SSTF decomposition of aP'^ involves surprising subtlety and concludes in Eqs. (|3|) and (Q) below. 

The boundary integral resulting from Eq. ( ^ ) is evaluated with (s', r', 0' , </>') coordinates and an expansion in terms 
of SSTF tensors. For simplicity we choose the boundary surrounding each source to be a surface of constant r' = tq; 
this boundary is spherical in a comoving frame of reference for non accelerating world lines and is appropriately 
Lorentz contracted in a frame wherein the world line is moving. Each boundary integral gives 

+ 0(G"+i). (37) 
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A Taylor series expansion of the Green function about the world line leads to the multipolar decomposition of the 
integrand in Eq. ( ^7|) on a two-sphere with constant retarded time, s'. But, the Green function is non-vanishing only 
on the past null cone from the field point Vap"'' picks up a contribution only where the past null cone intersects 
the three dimensional boundary, and this occurs for differing values of s. The point of intersection closest to the 
field point will have null-coordinate value s, but the point of intersection on the far side of the boundary will have 
null-coordinate value approximately s — 2rQ. Thus, the best two-sphere to use for the multipolar decomposition of 
VaP'^ is the one at s' = s — tq. And the Taylor series expansion for the Green function about the point z°'{s — s' + r') 
on a hypersurface which is orthogonal to the world line has x'"' = z°' + r'n"^ + 0{r'v°') with the right hand side 
evaluated at s' + r'; thus, 

G(a;" - x'") = J2 ^^^N'^VLGix"" - z^{s' + r')) 
/=o 

+ 0(G). (38) 
Now, with V^r'|s_ro — VJjr'j^ + 0{G) on the boundary and s' ^ s' — tq, 

\l+2 



(=0 ■ 



xVlG{x - z{s')) sine' d0' d(l)']ds' + 0{G''+^). (39) 

Below it is necessary that this integrand be in terms of SSTF tensors. To this end, is equal to a sum of terms 
each of which is a symmetrized outer product of projection operators /"i"^ and of the SSTF combinations N^; thus, 

[1/2] 
m— 

for a set of coefficients 6; m which are obtained in Appendix |c| and given in Eq. (C3). Now the substitution /°i°2 = 
-^0102 _|_ ya,iy0.2 g^jj^j ^j^g Qf (^^ allows part of the integrand of Eq. (|39|) to be rewritten for all x ^ Zs as 



[;/2] 

m=0 

x^lG{x-z,). (41) 



Application of the useful identity (B2) further transforms the u"^ . . . v""^"^ inside the integral into proper-time deriva- 
tives, and rearrangement of the summation ultimately yields 



oo „ 



;+2m+2 



47r(/ + 2m)! 



Lm=0 



ds'^" 

xWlG{x- z{s'))ds' + 0{G"+'^), (42) 



which has the ^ indices SSTF. 

Finally, the desired multipolar decomposition is 



oo „ 

VaK' = -11 K'J'{s)^lG{x - z,) ds 

+ 0(G"+i), (43) 



1=0 ■ 



where we define 



m— 

V av Elti{s - ra)N^ sine de d(j). (44) 



■' s,ro 
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The ^ indices are expUcitly SSTF, and K^^ has SSTF components P^, Q,^, Tl^ and defined from 

+ 7^f (.) + /^<'^'5^i>(.). (45) 
With this decomposition of V aP'^i hi hand, we re turn to the search for q"^ of the form given in Eq. (p4) which 



satisfies Eq. (p2|). This analysis closely follows §[II A. Let 



it = E / M:^\s)yLG{x - zs) ds (46) 

where the definition of M^^Hs) in terms oi A^. ■ -JCn is similar to Eq. (Isj). Now, when Vq?^^, as in Eq. (Hj), is 
added to V aP'^i^ as in Eqs. ( ^ ) and (jis]), and like terms arc matched up, the result of Eq. ( ^ is 

(47) 
(48) 
(49) 
(50) 

These form a set of coupled, linear, ordinary, inhomogeneous differential equations for A^^ ■ ■ ■ IC^ with sources involving 
7-"^ . . . 5^. And any solution to these equations gives a corresponding q^'' via Eq. (^4|) which, along with p'^^ determines 
5h^ and formally yields an improved, approximate solution to the Einstein equations. 

A particular solution to most of these equations results from setting A^, C^, V^, Tin, Jn and /C^ to zero, then 









' 71 


= 0(G"^ 


'1 


^ 




- ri n 




= 0(G"^ 


-1 








- 


= 0(G"^ 


-1 








~ <->n 


= 0(G"^ 


-1 



= -p-L =TZ ^ B^ 



gi:^Qt:, and £^ = 5^. (51) 

The general solution to Eqs. (|47|)-(|50|) is this particular solution plus any homogeneous solution for the A^ ■ ■ -l^n- 
And a homogeneous solution added in at the nth iteration is no different from starting the entire iterative process 
with a slightly different choice for the first order Ai . . . /Cf . 

The specific solution of Eqs. (^)-^(pO|) to be used should be determined by the physics of the interior. For example, 
the particular solution in Eq. ( |51| ) is easy to implement and leaves unchanged all of the mass and current moments, 
A^ and respectively. This may be loosely interpreted as the appropriate solution for a steady object and is not 
unreasonable as a choice for any astrophysically interesting source in a binary as long as tidal effects are unimportant. 
Also, this particular solution coupled with the simple choice for the first order moments of only the mass monopole, 
Ai, and current dipole, C^, being non-zero ought to be able to reproduce numerically the results of the higher order 
post-Newtonian analyses which are currently published and also contain no tidal effects. But, to find the solution 
appropriate for a tidally distorted star is more difficult — p"^ within the boundary creates tidal forces, distorts the star 
and changes all of the moments, A^ . . ./C^, in a manner which would need to be determined. Or, for a black hole 
within the boundary a pcrturbative analysis might be used to determine the appropriate solution of Eq. (p7[)-(p0[). 
In any event, a variety of different possibilities could be implemented; the actual choice made should be specific to 
the physics of the interior sources. 

But, the particular solution, above, fails for the low multipoles S„, JF^ and G'^ because these SSTF tensors don't 
exist and can't satisfy Eq. (|5T|). Thus, we still must contend with 

(52) 

A'.+S^+A'-n' = 0(G"+i), (53) 

(54) 
(55) 

These remaining ten ordinary linear differential equations have simple interpretations. An + T^n gives the rate of 
change of the mass monopole moment, and Vn is analogous to the rate energy flows into the source through the 
boundary. Cfj -I- jj^ is similar to the rate of change of spin angular momentum, and is analogous to the torque. 
B'^ + gives the rate of change of momentum of the source with respect to the world line, and TZ'^ is analogous to 
the force. A^ + 'Dn is closely related to the rate of change of the dipole moment caused by the momentum of the 
source with respect to the world line and by V^, which has no common Newtonian analog. 
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Eqs. js^ ) and ( |55| ) may be integrated as ordinary differential equations. Then, after a proper time s = 0{G~^), 
An, for example, will typically have grown large enough that the order of the approximation will have decreased by 
one. This is not particularly troublesome, and just implies that to obtain an approximate solution to the Einstein 
equations with E°-^{hn) — 0{G^) after a time s = 0{G~™') requires that n — N + m — 1. And for a binary system 
the approximation loses one order only when the fractional change in the mass of one of the components is 0(G"). 

But the physical interpretations of Eqs. (^3|) and (jsj) give cause for concern. So far in this formalism the world 
line has been given ahead of time. And the changing dipole moment and relative momentum of Eqs. (^) and (jsj) 
just reflect the fact that the true, physical source is moving with respect to the predetermined world line. But all of 
the moments of the source are calculated about the world line which does not necessarily follow the center of mass of 
the source. And as the source drifts away from the world line a rapidly growing number of multipole moments need 
to be monitored to adequately describe the source. This would be a disaster for any implementation. 

Thus at the n — 1 iteration, before getting to this stage, we should have made certain that the trajectory of the 
world line was chosen so that 

K-V'„ = OiG"+'). (56) 

Then with Bl = + 0(G"+i), Eqs. (||) and (||) are solved with all of A^, V\ and H\ being zero. And with no 
growing dipole moment the plethora of required moments is avoided. 

In the next section we show that Eq. ( pq) is essentially the (n — l)th order equation of motion of the source and 
generally necessitates an 0(G"~^) adjustment of the world line. 

At this point (if not previously) one might wonder what the effects of making a different choice for the radius of the 
boundary, tq, might be. For example, if tq were increased then Eq. ( |3^ ) shows that p'^^ would change by the addition 
of a homogeneous solution to the wave equation whose divergence would account for any consequent change in the 
Vn ■ ■ - Sn- Thus the change could be absorbed by q"^ . 



IV. EQUATIONS OF MOTION 



After the nth iteration of the field equations as described in §[IIB, both h'^ and z'^_i{s) are known. Before iterating 
the field equations again, it is first necessary to adjust the world line, z^_i(s) z'^{s), in order to enforce the nth 

order equation of motion, TZ^j^i — V^+i = 0(G""''^). 

Care must be taken to insure that this adjustment is accomplished while maintaining the order of accuracy of the 
current approximation to the metric. Thus, we require a satisfactory method for pulling the self-field of a source along 
a new world line. An aid in this task is the invariance of the Einstein equations under a Lorentz transformation, as 
formulated in §[IC. As we see below, the appropriate adjustment of a world line necessarily involves changing its 
acceleration but only by a small amount of 0(G"). And, pulling the self- field along a new world line is nearly, but 
not quite, accomplished by a Lorentz transformation. However, the retarded Poincare transformation, described in 
Appendix^ as a generalization of the Lorentz transformation, allows for a time-dependent boost and is still adequately 
behaved globally. The Einstein equations are not strictly inv ariant under a retarded Poincare transformation; but, as 
demonstrated for the scalar wave equation in Appendix D2, they are approximately so. And the retarded Poincare 
transformation is sufficient for the task. 

First in this section, we show how to implement the retarded Poincare transformation to pull the self-field of a 
source along a new world line in a manner that maintains the accuracy of the approximation to the Einstein equations. 
Then, we show just how the new world line is chosen to satisfy the nth order equation of motion. 



A. Adjusting 



A retarded Poincare transformation, described in Appendix adjusts a world line by defining a new coordinate 
system with 

y'^' =K''\{s,)x''+C'{s^), (57) 

where A"^ b{sx) is a matrix of the form of a Lorentz transformation and a function of the retarded time, Sx, at also 
f {sx) satisfies Eq. (p^). The world line in these new coordinates is 

zi =K'^\{s,_)zl_,+C'{s,), (58) 
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with the consequence, noted in Eq. (Dll), that 

vi ^ K'^\{s)v'^,^. (59) 

We show below how to choose A" b(s) to determine a new world line along which the equation of motion is satisfied; 
this has A° b(s = 0) = (5J^ , to match smoothly to the initial data, and A = 0(G") to keep the adjustment small. 
The field h"!^ is separated into a self field, h°^, and a background field, , 

K" = hf + /ilf (60) 



where hf^ contains at least the O(G^) part of the source whose world line is under consideration, and h'^ contains at 
least the O(G^) part of all other sources. The distribution of the remaining 0{G™;m > 2) parts of h'^ between h'^ 
and h"^ is immaterial. 

There is no unique way to pull the self-field h'^ along with the new world line. But if A" were constant then the 
usual Lorentz boost would be required. Thus, a natural choice for the self-field associated with a world line adjusted 
via a retarded Poincare transformation is 

hii^iy) = A^\A'' .h'^^ix), (61) 

and this seems particularly reasonable when A'' c ~ 0{G") and is small. Thus the new field, h'^^^{y), in the new 
coordinates is chosen to be 

htiiy) ^ A'^\A\ih^^{x) + <5f J^'/ig'(2/) 
= S^SMy) + [A'^'^A"' ah'^^ix) - <5f <5^'/i$f(y)] (62) 

where x is the function of y consistent with the inverse of Eq. (p7|). The background field h'j^ is the same function 
of ?/" as it was of while the self field is simultaneously pulled along and boosted by the time dependent Lorentz 



transformation, A° c{s)- This choice is consistent with the derivation of Eq. (75) below. 

But, we must show that the nonlinearity of E°-^{h) combines with A° to change Ef-^ only at 0(G"+^). The part of 
Kicwiy) ill square brackets in Eq. (^) vanishes when s is zero (because A"-' b{s — 0) — S^'), is proportional to h'^ and 
is, therefore, small and ^ sAHa — 0(sG"+^). Now, (^^wCy)) can be expanded about its value at 5^ h'^iv) 
and broken up into the parts 

E''''''{hii{y))=E'^''''[5:'6'^hl\y)] 

+ EfZ,,[A-' y ah'^^{x) - 5«'<5g'<(y)] 

+ -I-0(s2g2"+2), (63) 

where ii'Sncar denotes the linear part of the operator E°- ^ , from Eq. (|lj) ; and E'^ ^ is the part derived from r"'' which 
is still linear in A°'cA'''d/i5f (x) — S^' 5^ h'Aiy) but also depends upon S'^' S'^ h'^ (y) . The 0(s^G^"+^) terms remaining 
are at least quadratic in A° cA^ dh'^ix) — 5" ^Aiv)- observe how well h'^^^{y) satisfies the Einstein equations, 
we analyze Eq. (|6^) term by term. 

The first term is 0(G"+^) by assumption. 

The functional argument of the second term consists of the difference of two parts, each of which is 0(G), but 
whose difference is 0(sG"+^). The 0(G) piece of each part is a solution to the linear Einstein equations. And the 
difference of the O(G^) pieces of the two parts is only 0(sG"+^). Thus this second term is 0(sG"+'^). 

The third term consists of a sum of terms each of which is of the order of at least the product of 5^ h'^iv) with 
A''' cA''' dh^/ix) - 5"^ 5^dhA{y)\ the former is 0(G), the latter is 0(sG"+i). Thus the third term is 0(sG"+2). 

All together then 

E^'''[Kl{y)] = 0(G"+i) + 0(,sG"+2). (64) 
And we see that while s = 0(G~ ^) the error is 0(G"+^), after that the order of the approximation decreases by 



one in a manner similar to §IIIB. This is not a severe limitation on applications of this method to binary systems, 
where the radius of the boundary can be chosen to be of the same order of magnitude as the separation between the 
components, R. For a total mass of the system, M, and a typical speed, V, GM/R « (V/c)'^ = 0{G). The order 
of the approximation decreases by one only when sc/R — 0{G^^), and this occurs when sV ~ R{c/V)^. Thus the 
binary must orbit on the order of {c/V)^ times before the order of approximation decreases. 
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B. Adjusting the world line 



At this stage of an iterative step, we know z°_]^(s), h'^^{x) and h'^{x). This section shows how a change in the 

acceleration of the world line, effected by a retarded Poincare transformation with non- vanishing A" ^(s), insures that 
the nth order equation of motion is satisfied. 

For simplicity we assume that out of the first order moments, Ai . . . /Cf , only Ai is 0{G), and all of the others are 
smaller and O(G^). This is a reasonable physical assumption for astrophysical objects, and it should be clear how to 
generalize this analysis if one desires to examine, say, spin-orbit coupling by including terms here, as well. 

When the acceleration of the world line is changed, the dominant effect on TZ^-^ and arises from their parts which 
are linear in ^i, depend upon v°-, or and are O(G^). We call these parts 7?.^ and V^; and the next paragraph 
evaluates 7?.^ — Vj, with a conclusion in Eq. (|7|) . 

Both n\ and "Pji depend upon E^, the part of Ef_i which is linear in Ai and 0{G^). From Eqs. (|l2|), @) and 
it follows that 



(65) 



The evaluation of 7?.^ and V\ from Eqs. (Q) and ( p5|) requires the two integrals 



-5- / ^/avE^^rl smed0d(l) = -Ai{v'' + rov'') + 0{G^) 



(66) 



and 



(f \7arE^n'^r^ sine dOdcj}^ ^Airlv^v'^ + 0{G^), 



where use is made of the fact that VbV^ — —Vbv^ — 0{G'^). Now Eq. (|66|), along with h2r. 
(01), yields 



(67) 

l/(2m+ 1) from Eq. 



m=Q 

+ 0{G^). 



^2m 



J2m 



(2m + 1)! ds2r 



-iv' + rov'') 



(68) 



That the right hand side here is evaluated at s — ro is a complicating consequence of Eq. (44). Also, from Eq. (|6 
along with fei+2m,?ri = 3/ (2m + 3) from Eq. (C5), it follows that 



m=0 



r, 



2m+2 




(2m + 3)(2m+l)!ds2 



0{G^). 



These two equations combine to yield 



(69) 



(7^^ 



-{Aiv'')s 



2m+2 ,2m 

^0 a 



(2m + 3)! ds2 



r, 



2m+l J2m 


(2m + 1)! ds2m 



^2m+2 

'o 



(2m + 3) (2m + 1)! ds^' 



0{G^)- 



(70) 



the first two terms inside the summation come from the and ifi parts of 7?.^, respectively; the third term comes 
from V^. The two terms add directly, and the entire expression simplifies remarkably to 



{Tl\-V\)s^-A,Y. 



m jm 



0{G^). 



m! c?s'" 

Finally, the right hand side is a Taylor series expansion of —Aiv^ about s — ro but evaluated at s so that 



(71) 
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A 



(72) 



and both sides of this equation are now evaluated at the same s. The surprising simphcity of this last result might 
imply that a substantially more straightforward derivation could be found. 

A retarded Poincare transformation, with Eq. (p9|), effects the acceleration of the world line by 



(73) 



Now, we assume that A** ^(s) has been determined consistently with the equation of motion for all s up to some value 
So- Then both T^^+i ^'^^ ^n+i can be found from Eqs. (44) and (^ ) along with the temporary assumption that 
A''' a (so) = 0. And we define 



(7^: 



6' 

ri+1 



^n+l)so 



with A^'a(so)=0. 



(74) 



But, a value for A*''a(so) of 0(6"') changes (7^f,'+l 
Thus if we choose at sq that 



7^n+i)«o by -AA'''a<_i + 0(G"+2)^ from Eqs. (^ and (|7 



^iA\(so)<_i ^F" 



(75) 



then the adjusted world line will satisfy the nth order equation of motion at sq as welL 

This differential equation for A^ a is consistent with A^ q = 0(G"), as promised in §[VA, and gives three equations 

for A** a (the ^ index is orthogonal to w'' ); the remainder of A** a is determined by the requirement of Fermi- Walker 

transport Eq. (D21). 

With Eq. ([73D the differential equation can be rewritten as 



(76) 



which has the expected form for an iteration of the equation of motion with a residual force remaining on world 
line z„_i(s). It is not difhcult to show that at the first order this equation of motion is equivalent to the usual 
post-Newtonian result as presented by Bel et al. jl^ . 



V. BEHAVIOR AT FUTURE NULL INFINITY 



Now we reconsider the iterative procedure outlined in § [II with particular attention given to the limit of large r 



while s = t — r is held constant. Thus we consider the approach to future null infinity and show how to insure that 
the outgoing radiation propagates along flat space null cones which match up asymptotically with the null cones of 
the true, physical space-time. At every iteration an 0(G") gauge transformation, dX"^, and a small contribution, 
7"5j_]^ — 0(G"+^), insure that at large r, h"^ can be written as an expansion in inverse powers of r, times functions 
of retarded time, s, and angle, n, — in particular there are to be no Inr terms in this expansion. We refer to such an 
expansion as a proper expansion in inverse powers of r. 

We continue to use outgoing-null spherical coordinates, (s, r, 9, (f>), but now they are tied to a non-accelerating world 
line near the center of the binary system. 

At large r, h f^ admits a general multipolar decomposition just like that presented in Eqs. (^) and (|2^) and 
satisfying Eqs. (|2|)-(||) with the 0{G^) terms also being 0{r-^). Thus, 



hT^r-'xris,n) + Oir-') (77) 

defines Xi'', the dominant part of /i?'' at large r. But at second order, a difficulty immediately arises in evaluating pf'. 
Namely, Ef = 0{r-^) and Eq. (|3|) gives a r'^ Inr term to pf (Theorem 7.2 of Blanchet and Damour [|lo|). Such 
logarithmic behavior is the signature of a mismatch between the null cones of the background Minkowskii space and 



of space-time. Blanchet 12| shows that when the 0{r~^) part of E^^ is of a particular form, Eq. (85) below, then 
the logarithmic terms and the mismatch can be removed by a gauge transformation. Our analysis follows Blanchet 
1^ closely, except that we differ on a choice of gauge for hf' and that his analysis involves a clean separation of the 
powers of G, while our 0(G") terms contain further functional dependence on G. 

Analysis of the definition of r"^(/i) reveals that if kaxf' were zero then the offending part of E^'' would be easy 
to evaluate. Generally, kaxf' is not zero. In fact, it is straightforward, but not simple, to see that the restrictions 
(||)-(|8l) imply that 
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kaxf = -v\Ai+'Di)-{Bl + n1), (78) 

where Ai, B\, Vi and H\ now refer to the multipole moments of hf' as measured with respect to the non- 
accelerating center of the Minkowskii background geometry. With foresight, a Lorentz transformation removes the 
three-momentum, + ?ii)/4; a gauge transformation with V^A" — (this preserves the general form of hf'' and is 
discussed by Thorne |^) sets Vi = 0; and a second, preemptive, gauge transformation with A*^ (5-^1 Inr, —jAin°-) 
finally results in an hf'' whose lowest non-radiative multipoles take the simple form at large r 

hf « ir-Mlrfc^ (79) 

(Interestingly, this is the exact, at all orders, solution for /i'^'' for a Schwarzschild black hole in outgoing Eddington- 
Finkelstein coordinates [0.) With these gauge choices, ft-"'' is now in a form such that 

fcaxf = 0. (80) 

With this last result, the asymptotic behavior of Ef'' is dominated by r°''(/ii) and is of the form 

Ef = -r~^k''k^^2{s, n) -f Ef (81) 

where Ef = 0{G^r-^) and 

^2{s, n) = \xtxiab - IxlaXu = 0{G^). (82) 

We note that may be interpreted as being proportional to the effective energy density of the outgoing gravitational 
waves. 

Generally, this behavior for E^ occurs at every order. Rather than continuing a discussion with a focus on the 
second iterative order equations, we switch to the consideration of the nth iterative order and continue following 
Blanchet's ||l^ analysis closely. 

We iteratively assume that hf^^^i has a proper expansion in inverse powers of r times functions of s and n, 

h'i^_,^r~\':l,{s,n) + 0{Gr-^), (83) 
with 5xt-i = Xt-i - xt-2 = 0(G"-i), 

kaXn-i = 0, (84) 

and that 

= -r-^k''k''^r.{s, n) -f (85) 

where E^\ = 0(G"r-3) and = 0(G"). First we seek /if = hf_^+Shf^'' such that E''''{hn-i+Shn) = 0(G"+V-3) 
where this Sh^f differs from that of §111 by a gauge transformation and a small addition and, also, contains no Inr 
term at any order. Then we must reconsider the analysis of Eq. ( [l8| ) to account for the difficulties caused by the 
0(G"r-2) behavior of E^jLi- 

The initial task is to guarantee that E^^'^^^^^.{^hn) -Ef_^+0{C+^r-^) so that after the nth step E^'' will satisfy 
the incremented version of Eq. (|5|). The substitution of the Bianchi identity ( |l^ into Eq. (|35| ) results in 

VaK'=^ / K[E'tM^-^')]d' 

- ^ 1 E^i, ivc'ri^ ~ rl) G{x - x') d\'. (86) 

The second integrand is 0{G'^'^^r~^), thus the integral is 0(G""'"^) and has a proper expansion in inverse powers of 
r the leading term of which matches an outgoing solution of the homogeneous vector wave equation. In §|l|, q"^ was 
chosen to cancel just the 0{G") part of VaPn'''^ now, qf in Eq. ( ^ ) can be chosen at no additional expense to cancel 
this 0{G^~^^r~^) leading term from the second integral as well. This results in 

Va(K' + Qn) = r-'C+ii^, n) + 0(G"+V-3) (87) 

for some vector ('^_^_-^{s , n) = 0(G"+^). Further, if we choose 7°'^^ = 0(G"+^r^^) such that 
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fca7,'?ti=^"'C.+i+0(r-3), 
Then + q'^ + 7^^^ is nearly what we seek for Sh"^ . With the definition 

Tjab — ab , ab , ah I oc\\ 



it follows that 



Also 



Va-fff = 0(G"+V-3). (90) 



= 0(G'"+V-3), (91) 



so that 



V^ijf = i;,^^, + V'7,'J^i = + 0(G"+V-3). (92) 



Both Eqs. (|0|) and (||) are used below. And E^^^^^{Hn) = -^^,'^'^1 + 0(G"+V-3) as required, so iJ,f would be the 
choice for 5h'^ except for the Inr behavior of p"^ . 

linear 



An 0(G") gauge transformation resolves this difficulty and leaves E^^^^^ unchanged. Blanchet Q shows with his 



Lemma 2.1 that if 



then 



A" 



'^n{t,n')dtG{x~x')d^x' (93) 



and 



dK' = ^J ^^ni-s', n') Gix - x') d'x' 

+ 0(G"r"i) (94) 

VaSAf ^-TT^r r *„(s',n)ds'. (95) 

Thus, from Eqs. (^3|), ( p5| ) and ( p^ the combination pj^'' + 9AJ^^ is a proper expansion in inverse powers of r. And, 
with 

oK =Pn +Qn +0Xn + ln+1 (96) 

and = ^n-i + ^f^n''' it follows that Sh'^'' has a proper expansion in inverse powers of r, and Eq. (^3|) holds with 
n — 1 —> n. 

The iterated versions of Eq. (E4|) and (|5|) remain to be checked. From Eqs. (po| ) and ( |95| ) it follows that V aSh'^i — 
0(G"r-2). But with = r-^xf (s, n) + 0(G"r-2), it must also be that 

VaSh'i^ = -r-ifc.Jx:^'' + 0(G"r-2). (97) 

Hence kaSxl^ = 0; and, for n > 1, Sx'!'!' is zero on the initial hypersurface, where hf' matches smoothly onto the 
initial data, so kaSx^ = always, and Eq. (84) holds with n — 1 — > n. 

Finally, the analysis of Eq. ( [l8| ) is modified by the presence of dX'^'' and "fn+iJ but, 

E-\h,,) = [E-\K,) - E-\K^i)] + E''\K-i) 

= V'^VcK" + V^V.if™ - v'^^'^c^dH^f - IQA-r'^^'iK) - r'^''(/i„-i)] + 0(G"+V-3) (98) 

where the second equality follows from the definition of E°'^(h), Eq. (|l^). And the third equality is a consequence 
both of 0(G") gauge invariance of i^iJnoar ^^'^ ^1- (92). Now, Eq. ( pO[ ) and the application of the iterated 

versions of Eqs. (||) and {pi) to the definition of r°''(/i„) imply that E''^{K) is of the form of Eq. (||) with 
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*„+l(s,n) = ^ixtXnah - Xn-lXn-l ab) 

~ jiXnaXnb ~ Xn-laXn-lb)- {99} 

An iterative step is thus formulated in a manner which leaves h'^'' expressible as a proper expansion in inverse powers 
of r with Eqs. (|8^)~(p5|) holding at every iteration. And the gravitational waves asymptotically expand out along 
constant s surfaces, so the outgoing null cones of the true space-time metric asymptotically match up with the flat 
space null cones. 



VI. CONCLUSIONS 



We have given a prescription for iteratively improving an approximate solution to the Einstein equations which 
could be carried out, by computer, to any order. The lowest order approximation is just the familiar linearized 
approximation of general relativity. 

The description of the iterative process in the text was intended to follow a logical order to provide motivation for 
each part, in turn, of one full iterative step. However, in practice the chronological order is slightly different. Thus, 
we now summarize the entire process with a brief chronological description of the procedure. 

For the initial step we choose ft."'' to be the algebraic sum of two terms like Aiv°'v^/r, one for each of the two 
sources. Let this hf' be considered h^:^ and let the world lines be z^'Li- 

In order to iterate the equations of motion, first find E'^J' on each boundary and then find TZn+i ^nd V^^^i from 
Eqs. ( PI ) and (^. Now, Eqs. and (D21) determine the A"^';, which adjusts the world fine to and modifies h'^'' 
to h^i^ with Eq. (H) while preserving E^n^^) = 0(G"+i). 

Next the field equations are iterated by using Eq. (|33|) to determine pj^+i and Eq. (|4^ ) for gj^'^^, with the [n + l)th 
order moments satisfying both Eqs. (^7|)-(pO|) and also appropriate conditions determined by the physics of the 
sources within the boundaries. The determination of p'l^_^_i involves solving the inhomogeneous wave equation for all 
independent components of the symmetric tensor. This is the single, computationally intensive part of every iterative 
step. 

At this point E'^^ihn+Pn+l + gn+i) = 0(G"+2). But to preserve proper behavior at future null infinity with the 
outgoing null cones of flat space-time matching up asymptotically with those of the true, physical space-time, a gauge 
transformation, dX^j^^, from Eq. ( p^ ) is needed. 

In preparation for the next iteration g^+i should also cancel the 0(G'"+^r^^) contribution from the second integral 
in Eq. ( ^6|) . And h^_^_i should be changed by a small correction, ~ 

0{G"+^r-^), which satisfies Eq. (||). Now, 
/i^^i = ft"j'' + Pn\i + + dX^^\i + 7^5+2 completes one fuU step of the iterative procedure. 

The freedom of this iterative process from the restriction of the harmonic gauge may provide an important aid to 
its implementation. For example the first order /i"'' might be chosen to be the sum of two terms like Ak'^k^ /2r. Then 
individually each term would be the exact Schwarzschild geometry, if the source were not accelerating. And E°'''{hi) 
would consist only of linear terms dependent upon the acceleration and cross terms between the two sources. This 
choice for hi'' would already be an accurate approximation to two Schwarzschild black holes even near the past event 
horizon of one hole where E°'^{ hi) ^ M/R, with M and R being the mass of and distance to the companion hole. 

One weakness described in § IIIB| stems from the inability to treat the conditions at the inner boundaries in a 
straightforward manner. In problems whose focus is on the emission of gravitational radiation from binary systems, 
this is a difficulty only when the system is tight enough that tidal deformations are important. To include tidal effects 
of any sort it is necessary to solve the internal problem, as described for example by Damour , and thus to obtain 
the specific solution of Eqs. (^7|)-(^0|) which matches the physics of the problem. 

Work in progress applies similar methods to a Schwarzschild background geometry. In this case it appears that 
perturbation analysis ought to yield boundary conditions which can be properly imposed at the event horizon. This 
extension will provide a better method for the analysis of black hole binary systems. 
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APPENDIX A: THE RETARDED GREEN FUNCTION. 



The retarded Green function, 

G{x~x')^20{x° -x'°)S{n), (Al) 
where 57 is the square of the flat-space interval between two points, 

n{x,x') = rjabix" - x"')ix'' - x"'), (A2) 

is a solution of 

W''WaG{x-x') = -4TTS'^{x-x'). (A3) 
For a generic problem we wish to find a particular solution of 

V^/i"^ = -p"^ (A4) 
where p""^ is a multipolar skeleton source on z°(s). Thus, if 

p^^x) ^Att j M"^^(s)VL<5''(a; - z,) ds, (A5) 

then 

= I G{x~ x')Vl I M''''^{s)5\x' - zs) ds d'^x' , (A6) 



with V^, being the derivative operator with respect to x''^ . 

After integrating by parts I times, changing the derivatives to be with respect to x"" and integrating over a;'°, we 
have 

h''\x) = j hr^^{s)VLG{x- z,)ds, (A7) 

or after withdrawing from the integral we have 

/^'•^(x) = Vi[r-iAr^^(s,)]. (A8) 

While this result may appear quite familiar, it is important to remember that it holds for a source which is moving 
along some accelerating world line and has time-changing multipole moments all the while. The consequent radiation 
results from both the acceleration as well as the varying multipole moments. 



APPENDIX B: A USEFUL IDENTITY 

A useful identity is 

j f{s)v''VaF{x ~Zs)ds^-J f{s)^/{x - Zs) ds 



= -[f(s)Fix-z,)\°^^ + J ii^Fix-Zs)ds. (Bl) 

For our applications this integral is over all proper time, s, and F{x — Zs) involves the retarded Green function and is 
zero except where the past null cone from a;" intersects the world line. With these conditions the contribution from 
the limits of integration is always zero, and 

f{s)v-\/aF{x - Zs) ds= I ^^F{x - Zs) ds . (B2) 



This identity is particularly useful in the reduction leading to Eq. 
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APPENDIX C: SSTF DECOMPOSITION OF A^^ 



By repeated subtraction of the trace parts, iV^ is expressed as 



[1/1] 

m=0 



for some set of coefficients 6/.m- Contraction with fai^iai and use of Eq. ( |Cl| ) with I —>■ I — 2 yields 

{2m + 2) {21 - 2m - 1) 



l{l-l) 

It is clear that bifi = 1, and with elementary methods we find that 



-&i,m+l- (C2) 



^ /!(2?-4m+l)!! 

2™m!(Z-2m)!(2Z-2m+l)!!' ^ ' 

Some special values which are of use for determining TZ^ and are 

62™,™ = l/(2m+l) (C4) 

and 

= 3/(2m + 3). (C5) 

APPENDIX D: RETARDED POINCARE TRANSFORMATIONS 

The retarded Poincare transformation is a little known method for relating outgoing-null coordinates associated 
with different world lines. This transformation is a mapping from one flat space-time to another, which transforms 
one world line, z°'{s), into a second while preserving the values of the scalar fields s and r; also, the future null cone 
of each event on the first world line is mapped to the future null cone of the corresponding event on the second world 
line. Strictly speaking this transformation is a diffeomorphism from the causal future of z°'{s) (that is from the set 
of all events which can be reached from z"(s) by a future directed, non-spacelike curve) onto the causal future of the 
second world line. This technicality is required to allow for the possibility that one world line has constant acceleration 
in the distant past, and its causal future is, thus, not the entire Minkowskii spacetime. 

We first describe the mathematical form alism of the retarded Poincare transformation and then give an application 



which is closely related to the analysis of §IV A. 



1. Mathematical formalism 

We start with a given world line z°(s) in Minkowskii space covered with the usual Minkowskii coordinates, x", and 
define a coordinate transformation by 

y-' ^A-\{s.,)x'+e'{s^), (Dl) 



where A" b and are functions of s^, explicitly, and of x"", implicitly, and i^" (s) satisfies Eq. (D6), below. The 
matrix A" is a time dependent Lorentz transformation, i.e. it is a matrix of the general form of a Lorentz boost and 
a rotation, as described by Misner et al. JlGf , but with the boost and rotation parameters being functions of s; also, 
A^'c' is the matrix inverse of A° Thus transformation (Dl) reduces to a Lorentz transformation if A° f, is constant 



and = 0; thus, this transformation is a time dependent generalization of the Lorentz transformation which is 
reasonably well behaved in a global sense. 

In this section a prime on a base letter identifies a geometrical object which is most naturally discussed in the 
coordinate system; a prime on an index refers to the components of a geometrical object in the coordinate system. 

We also define 
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^a'h' ^ A°',A^'d7f (D2) 

From the algebraic properties of Lorentz transformations we know that the y°' components of 77'° ^ are (—1,1,1,1) 
on the diagonal and zero elsewhere. While the usual flat Minkowskii metric of the coordinate system, it 

is not the tensor equivalent of rf^ with the coordinate transformation (Dl), because 



'i.'-i.r + r,"'- (D3) 



We further define 



7?:,,, ^ k\,t.\,^,i. (D4) 



Then is the matrix inverse of t]'^ because of the usual algebraic properties of Lorentz transformations. And we 
raise and lower primed indices on primed tensors with Vj^ and v(y^i\ however, with two different metrics at hand we 
rarely raise or lower indices implicitly. From these definitions it follows that 

T(^,yk''\^^,,k\. (D5) 

along with some index variations of this equation. 

For the given world line, 2;''(s), we choose (s) so that 

t' - -k'^\z\s)- (D6) 

this uniquely determines (s) up to the addition of a constant vector. This choice is motivated below, after Eq. 

The transformation of the components of tensors is governed by 

^K'^'b-rk'^' Xh. (D7) 

The inverse transformation is 

h.\,-rM'd'k''kar, (D8) 



the derivation of this inverse involves some of the results d eriv ed below. 

We are free to consider the coordinate transformation (Dl) as a diffeomorphism from one manifold to a second. 
Then 77"^ is a flat metric on the x"" manifold, and ry'" ^ is a flat metric on the y° manifold. With this point of view, 
the world line z°'{s) is mapped to a world line on the manifold by 

z"'\s)=K''\z\s)^C'(s). (D9) 

Then 

i'^'is) = K''\z\s) + k'^\z\s) + r'; (DIO) 



and with Eq. (D6), the four- velocities are related by 

v'"' = = A°'fc^;^ (Dll) 

and it follows easily that jy^/^/w" = — 1 demonstrating that s is the proper time for the world line z''^ as well. 

The future null cone structure of the world line is preserved under the transformation in the sense that the future 
null cone of the event z°'{s) is mapped onto the future null cone of the event (s). For a proof consider the square 
of the interval between a generic point on the y"' manifold and a point on the world line, 

(/,z-')^(/-z-')(/-z"'')r,l,,, 
= n{x\ z^) + 2A'^\{x'' - z'W,,,, [C' is.) - r' (s.)] 

+r' is.) - r' is.)Wa'c' r' («.) - r' («.)], (D12) 
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where the second equahty follows from Eqs. (Dl) and ( D10| ). Thus, if a;° is on the future null cone of z"" then Sx = Sz 



and r2(a;°, z"") = so that ^' {y'^ , z'"' (s)) = 0; it follows, then, that j/" is on the future null cone of z"^ as determined 
by the Ty^/f,/ metric. 

We define a useful scalar field, r', similar to r in Eq. (^, by 



(y" - z"' ) = -r^.dv'lx'' - z'^isx)] = r. (D13) 



In other words, r'{y{x)) = r{x). 
Finally, from 



= iA'''a-rA'''ck'ka)K,, (D14) 

and contraction with k°' reveals that fc^'A'' afc^, = 0; with some effort it also follows that ak'i,,k°' — and that both 

ka^A^'ak',,, (D15) 

and 

k" ^ A-'b'k'^' . (D16) 

With these results at hand, we simplify and summarize the notation. The x"^ manifold and the manifold have 
some similar structures which are distinguished by a prime for the structure on the manifold. Examples are 
z'" and 77^/ J,/ which are similar to and rjab- But if we consider the mapping from to to be a coordinate 
transformation then wc can den ote the components of 77' in the coordinate system as 77^^, and as discussed above 



^afc 7^ ""lab- But from Eq. ( D14 ), it does follow that ka = k'^, so that ka and k'^, are the same vector in different 
coordinates; hence we leave the prime off the base letter k; and an index on k can be raised or lowered by either i] 
or 77'. Additionally, both s and r evaluate to the same scalar fields on the two different manifolds — and we leave the 
primes off these fields as well. 



Thus we see that the retarded Poincare transformation described in Eq. (Dl) is a diffcomorphism which maps one 
world line into another while preserving its future null cone and the values of the scalar fields s and r. Associated 
with each world line is a distinct flat metric, rj""^ or 77'" , which is of the usual Minkowskii diagonal form in the 
appropriate (resp. or 77° ) coordinate system. We find it most convenient to be able to move easily between these 
two manifolds. 

It is not difficult to show that the composition of two retarded Poincare transformations can be described as a 
single transformation, and also that any retarded Poincare transformation has a unique inverse. 

For any world line, z'" (s), there are many retarded Poincare transformations from z"(s) — (s, 0,0,0) to z'" (s) 
which also rotate the coordinate basis vectors, . But, if the are Fermi- Walker transported along z'"' (s), then, 
as we now show, the transformation is unique up to an initial rotation of the basis vectors. Fermi- Walker transport 
requires that 

ef =~il/^b'e^, (D17) 

where 

n""'^, = v""' v'^, - v'"' v'^, . (D18) 
In addition, for a retarded Poincare transformation the coordinate basis vectors must also obey 

ef'^A'^'bel (D19) 



The substitution of Eq. ( D19 ) into Eq. ( D17 ), along with the orthonormality of the basis vectors, results in 

A'^'b = -n"'\.A''b, (D20) 



which has a unique solution, given suitable initial conditions. This last equation, along with Eqs. (D6) and ( |DS| ), 
determines the retarded Poincare transformation that maps z"(s) = (s, 0, 0, 0) into z'" (s) with Fermi- Walker transport 
of the basis vectors. 

In the more general circumstance that z"(s) as well as z'" (s) are arbitrary world lines, Fermi- Walker transport 
requires that 

A'^'b = -^Q.'"' c^A-^'b + A'^'.r!^,, (D21) 



where VL'^b is defined as in Eq. (D18|) but with the primes removed. 
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2. An application 



An interesting application of the retarded Poincare transformation, related to the analysis of §[VA, involves a 
scalar field which satisfies the wave equation with a 2'-pole source moving along some given world line and with no 
incoming radiation at infinity. 

For a time dependent source at rest in the x manifold, a simple expression for the scalar field is 

^ (-!)'(/ + fc)! ii'^)M^{s)NL 



k=0 



where A^^(s) are the retarded-time dependent SSTF 2'-pole moments of the source, and the prefix superscript of 
M'" denotes differentiation with respect to s; this equation is given by Thorne Q in Eq. (2.53a). 

More generally, let a retarded Poincare transformation generate a different world line, z'" (s), for the source on 



the y° manifold via Eqs. (Dl) and (D9); then the solution of the scalar wave equation, retaining the same 2 -pole 



moments as measured by a nearby comoving observer, is easily written as 

V''(y)=V'i,(A^^'(s)/r), (D23) 

where r is defined as in Eq. (||), and A4^ = k-M^ . But for an accelera ting world line the evaluation of this right 
hand side for tfj' is much more complicated than the right hand side of Eq. ( |D22| ): it involves some terms with up to I 
derivatives of the velocity and others containing u'; and, even though the sour ce is a 2'-pole, ijj' has moments for the 



radiation all the way from a monopole up to a 2^'-pole. Thus, although Eq. ( D23 ) looks simple, it is actually quite 
difficult to evaluate for an accelerating world line. 



Now, consider a new scalar field on the y° manifold defined in terms of the scalar field of Eq. (D22) on the x 
manifold by 

^ ^{x); (D24) 
this is not a solution of the scalar wave equation on the manifold, but for small accelerations it is nearly one. In 



fact, inductive evaluation of Eq. (D23) along with Eqs.(W) and (g) and the invariance of s and r under the retarded 



Poincare transformation, shows that 



<ow(2/)-^'(2/)(l + 0(A)), (D25) 



and this is uniformly valid even for large r. 

A more general analysis reveals that in a similar manner a retarded Poincare transformation can take a scalar wave 
solution for a multipole source moving along one accelerating world line and generate an approximate solution for the 
same multipole source moving along a different world line, as long as the difference in the accelerations of the world 
lines is small. Note that this final limitation restricts neither the total acceleration of the trajectory nor the size of 
the boost allowed in changing the world line, as long as the proper time derivative of the boost is small. In this same 
manner Eq. (^TJ) gives a new self field, h°^-acwi when the world line of the source is changed via a retarded Poincare 
transformation. 
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